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ABSTRACT 

We report the discovery of near-sinusoidal radial velocity variations of the GOV star p CrB, 
with period 39.6 days and amplitude 67 m s _1 . These variations are consistent with the existence 
of an orbital companion in a circular orbit. Adopting a mass of 1.0 Mq for the primary, the 
companion has minimum mass about 1.1 Jupiter masses, and orbital radius about 0.23 AU. 
Such an orbital radius is too large for tidal circularization of an initially eccentric orbit during 
the lifetime of the star, and hence we suggest that the low eccentricity is primordial, as would 
be expected for a planet formed in a dissipative circumstellar disk. 



Subject headings: planetary systems — stars:low-mass, brown dwarfs — stars: individual 
(p CrB) — techniques: radial velocities 



1. INTRODUCTION 

Within the past 18 months, our knowledge of low mass companions to stars has exploded, so that it 
is becoming possible to carry out comparative studies of their properties such as mass, orbital period, and 
orbital eccentricity. This information is opening the door to a better understanding of the formation and 
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evolution of planets and brown dwarfs, including the planets in our own solar system. However, the small 
number of objects with known orbital parameters does not yet permit creating a reliable picture of the 
origin and evolution of such objects. 

For this reason it is extremely important to increase the number of known low-mass stellar companions 
and also gain information on their orbital properties. We have developed the Advanced Fiber Optic Echcllc 
(AFOE) spectrograph in part for this purpose. Using this instrument we have begun an intensive program 
of monitoring about 100 solar-type stars brighter than V ~ 6.5 with a precision of order 10 m s _1 ; this is 
sufficient to detect the reflex orbital velocities of any of the substellar companions detected by precise radial 
velocity techniques to date. In this Letter we report the detection of periodic radial velocity variations in 
the star p Coronae Borealis, which we interpret as due to a Jupiter-mass planet in a near-circular orbit 
about the star. 



2. RHO CORONAE BOREALIS AN OLD SOLAR-TYPE STAR 

In order to understand the properties of an orbiting companion to p CrB (HD 143761, HR 5968), it is 
necessary to estimate its mass, effective temperature, and age. The evidence argues for p CrB being an old, 
solar type star, with mass close to that of the Sun. 

The spectral type of p CrB is reported to be cither GOV (Jaschck, Conde, & de Sierra 1964) or G2V 
(Roman 1955, Roman 1952). Reported B — V color indices average about 0.61 (Roman 1955, Naur 1955, 
Argue 1963), i.e. slightly bluer than the Sun, for which reported values of B — V range from 0.63 to 0.66 
(e.g. Gray 1992, Cayrel de Strobel 1996, Taylor 1994). The effective temperature T cff is about 5760K as 
determined from the infrared flux method (Gratton, Carretta, & Castelli 1996 and Alonso, Arribas, & 
Martinez-Roger 1996), although Gray (1994) concludes from line-depth ratios that T c ff = 5868 K. 

The luminosity of p CrB is larger than the Sun's, suggesting that p CrB is more evolved and thus 
older. This is inferred as follows. Its V band apparent magnitude is close to V — 5.41 (Rufener 1976, Naur 
1955, Argue 1963, and Roman 1955). A recent (pre-Hipparcos) estimate of its parallax (Gliese & Jahreiss 
1991) is 7r = 60±6 mas, giving a distance of 16.7±1.7 pc. These values imply an absolute visual magnitude 
My = 4.3. Using a Sun-like bolometric correction of B.C. = —0.07, this gives a luminosity of L = 1.61 L Q . 
Adopting an effective temperature equal to that of the Sun (5777 K, Cayrel de Strobel 1996), we obtain a 
radius of R = 1.27 R Q . 

Based on evolutionary H-R diagrams of Bressan et al. (1993), the color and luminosity of p CrB are 
consistent with a 1.0 solar mass star about 10 Gyr old. The surface gravity was determined by Kiinzl et al. 
(1997) and Gratton et al. (1996) to be log(g) = 4.23 and 4.11 respectively, which is also consistent within 
the uncertainties of the fitting procedure with an evolved solar mass star. 

Henry (1997) reports that four years of precise photometric data for p CrB show no significant 
periodicities between 1 and 50 days, and constant seasonal means to within 0.00017 magnitudes. The 
implied absence of spot activity is consistent with an old age for p CrB. 

The relative chromospheric inactivity of p CrB also supports its high age. Its Ca II chromospheric flux 
is below that of the Sun and shows no significant seasonal variations (Baliunas et al. 1995). In spite of the 
smaller Ca II flux, the rotation period of p CrB as predicted from that flux is about 20 days (Noyes et al. 
1984, Sodcrblom 1985), which is slightly shorter than the 25.4 day sidereal rotation period of the Sun. The 
faster predicted rotation results from the fact that p CrB has slightly bluer B — V color index than the Sun. 
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There is some observational evidence for a rotational modulation of the Ca II flux with about a 17-day 
period (Baliunas et al. 1996), although that evidence is marginal (Soon 1997). 

Further evidence that p CrB is older and more evolved than the Sun is provided by the metal deficiency; 
reported values of [Fc/H] average to —0.19 (Kunzl et al. 1997, Hearnshaw 1974, Wallerstein 1962, Alexander 
1967). The ultraviolet excess S(U — B) = 0.06 (Hearnshaw 1974) also indicates a slightly less than solar 
abundance. In addition, p CrB is a high proper motion star with a motion out of the galactic plane of 28 
km s _1 (Cayrel de Strobel 1996), indicating that it may be a member of the old disk population. 

The above information, taken all together, suggests that p CrB is probably of near solar mass but 
older, near the end of its main sequence lifetime or possibly already starting hydrogen shell burning (cf. 
Cayrel de Strobel 1996, figure 14). For the purposes of this Letter, we shall adopt a nominal mass for p CrB 

Of 1.0 Mr,,. 



3. OBSERVATIONS 

3.1. Instrumentation and Reduction Procedure 

Data on p CrB were taken with the Advanced Fiber-Optic Echelle (AFOE), located at the 1.5-m 
telescope of the Whipple Observatory. This instrument is a cross-dispersed, fiber-fed echelle spectrograph 
designed specifically to perform precise stellar radial velocity measurements (Brown et al. 1994); its 
construction and operation are a joint project of the Smithsonian Astrophysical Observatory, the High 
Altitude Observatory, and Pennsylvania State University. The spectrograph covers a wavelength range of 
about 150 nm between wavelengths of 392 nm and 665 nm, at a typical resolution R = X/5X of 50000. 
The combination of a fiber feed, mechanical stability, and careful thermal control minimizes drifts in the 
velocity zero point and also results in a very stable instrumental resolution profile. For these observations, 
wavelength calibration was provided by an absorption cell filled with iodine vapor (Horner 1996), which 
superposes a dense spectrum of molecular iodine (I2) lines upon the stellar spectrum in the wavelength 
range between 500 nm and 620 nm. When operated in this mode, the spectrograph can achieve Doppler 
precision of better than 10 m s _1 on a 5.0-magnitude G star in an observation time of 10 minutes (Noyes et 
al. 1997). 

The method of data reduction is generally similar to that described by Butler et al. (1996), although it 
differs in details. 1-D spectra are extracted from the 2-D echelle images, corrected for scattered light, and 
divided by the flat-field spectrum of a continuum lamp. For each of six spectral orders containing strong 
I2 lines the stellar radial velocity is determined independently, by modeling the observed star-plus-iodine 
spectrum as the product of a Doppler-shifted, high signal-to-noise-ratio reference spectrum of the star alone 
and a very high spectral resolution spectrum of the I2 absorption cell in the laboratory reference frame; the 
latter was previously measured at the McMath-Pierce FTS spectrograph at the National Solar Observatory. 
The model incorporates the sought-for relative wavelength shift between the spectrum of the star and that 
of the iodine cell, detector shifts due to instabilities within the spectrograph, and a parametrized wavelength 
solution and instrumental resolution profile. The instrumental profile shape is allowed to vary slowly along 
the order. 

All of the parameters describing this model are adjusted so as to minimize the rms difference between 
the observed and modeled spectrum. The six independently deduced Doppler velocity shifts for each order 
are then combined with equal weights to provide the final velocity shift for that exposure, while the internal 
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scatter amongst the different orders gives an estimate of its uncertainty. Each observation normally consists 
of three consecutive exposures, in order to correct for cosmic ray contamination; typical exposures have 
continuum signal-to- noise ratio per pixel of order 100 to 150. The results of the velocity estimates for the 
three exposures are corrected for motion of the telescope relative to the solar system barycenter, and then 
combined to provide the final measure of the stellar radial velocity shift and an estimate of its uncertainty. 



3.2. Observations and Results 



Observations of p CrB were obtained on 41 nights between 1996 May and 1997 March. A tabulation 
of the data, showing times of observation, derived velocities, and estimated uncertainties, may be found 



on the World Wide Web at |f tp : //cf a-f tp .harvard. edu/pub/af oe/rhoCrB/| . The data, together with 
uncertainties, are plotted in Figure 1. 
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Fig. 1. — Top: Radial velocity of p CrB (solid circles), with three 2.5-cr rejected points shown as crosses. 
Dashed line is an orbital solution using the parameters listed in Table 1 . Bottom: residuals from the orbital 
solution. 

We fitted a Keplerian orbit to the data using a least-squares fitting process. The orbital parameters 
are given in Table 1. They result from an iterative 2.5-er rejection fitting procedure (where o is the rms 
of the residuals to the fit); three points were rejected on this basis, as shown by crosses in Figure 1. 
We note that the orbital parameters resulting from a fit without any rejection, or from a 3-cr rejection 
procedure which eliminated only one point, differ from the values listed in Table 1 by less than their formal 
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uncertainty. We verified that the elements we derived are the ones providing the best fit to the data by 
using a x 2 -nfinimization routine based on a genetic algorithm (Charbonneau 1995). This algorithm uses a 
form of directed Monte Carlo search that is specifically designed to locate the global minimum, regardless 
of where it lies in the parameter space. The fit to the next most likely period (at 76.8 days) has residuals 
with reduced % 2 of 3.3, corresponding to an extremely small probability that this is the correct period. 

The period and amplitude of the fit are well determined, but because of the small eccentricity, the 
longitude u> and epoch T of pcriastron have large uncertainties. However, the epoch Transit of passage of 
the companion across the line of sight to the star is reasonably well determined, as shown in the table. This 
value is included to facilitate possible searches for evidence of a transit across the disk of the star, should 
sin i be close to unity. However, Henry (1997) reports that his 208 photometric data points over the past 
four years, when phased at the period in Table 1, make a transit unlikely. 

The derived eccentricity (0.028 ± 0.040), while close to zero, has sufficient uncertainty that a value 
as large as 0.1 cannot be ruled out. As noted in Section 4 below, whether the eccentricity is significantly 
different from zero is an important question; more observations, especially with complete phase coverage, 
are needed to pinpoint the value of the eccentricity more accurately. 

The dashed line in Figure 1 shows the Keplerian orbital fit from Table 1. Residuals from the fit are 
shown in the lower part of Figure 1. Figure 2 shows the same data phased according to the orbital fit. The 
periodic variation of radial velocity is evident. 
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Fig. 2. — Radial velocity of p CrB (circles), phased according to the orbital elements listed in Table 1. Solid 
line is the orbital solution. 
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Part of the residuals from the orbital fit may be attributable to long-term instrumental drifts of 
uncertain origin. In Figure 1 , such drifts are suggested especially by the groupings of residuals for the four 
1996 observing runs (HJD < 2,450,450); however the smaller rms residuals (only 5 m s _1 ) of the data points 
for the 1997 January, February, and March observing runs indicate that the problem may have largely 
disappeared after an instrumental realignment in early January of 1997. Here we take the conservative 
approach of not correcting for long-term drifts in the orbital analysis. 



4. DISCUSSION 

The data plotted in Figure 1 clearly indicate a nearly sinusoidal variation in the measured radial 
velocity of p CrB, with period of 39.6 days and amplitude of 67 m s _1 . The obvious interpretation is that 
p CrB is orbited by a low mass companion, which produces a reflex motion of p CrB itself as they revolve 
about their common center of mass. If the mass of p CrB is 1.0 M Q (Section 2), the minimum mass m<2 sim 
of the companion is 1.1 Jovian masses, and the orbital radius a-i is 0.23 AU. 

Other interpretations could be considered. It is possible, although unlikely (Soderblom 1985), that in 
spite of the approximately 20 day rotation period implied by the mean Ca II flux level, the actual rotation 
period of p CrB is 40 days. In that case, an apparent periodic variation of radial velocity could in principle 
be caused by changes of the photocenter of the rotating star by periodic passage of a large spot across its 
disk. However, an amplitude as large as we observe would require about 3% spot coverage, which would 
produce a photometric modulation of order 0.03 magnitude. Henry (1997) reports an upper limit of 0.0002 
± 0.0001 magnitude for any variability at the period of 39.65 days. This result effectively rules out the 
above interpretation. A rotating magnetic region with a convective blueshift which somehow produces 
no photometric signature could be considered, but the more magnetically active Sun shows no spurious 
radial velocity variations as large as one-tenth the amplitude we observe on p CrB (see McMillan et al. 
1993). One might postulate a non-radial pulsation which does not produce any disk-integrated photometric 
modulation, as Gray (1997) has suggested to explain the 4.2-day periodicity observed in the radial velocity 
of 51 Pegasi (Mayor & Queloz 1995). However, the 40 day period of the radial- velocity variations in 
p CrB is extremely long compared to the dynamical timescales characterizing acoustic-gravity modes in 
Sun- like stars (roughly 1 hour). Moreover, the amplitude of 67 m s _1 is enormous compared to the 1 m s _1 
amplitude of the combined solar "5-minute" p mode pulsations, or of the (so far undetectably small) solar g 
modes. Furthermore, the abovementioned stringent limits on photometric variation are extremely difficult 
to reconcile with radial or low-degree nonradial pulsations at an amplitude approaching the observed 
velocity amplitude. For these reasons, we take the view that the radial velocity variation observed in p CrB 
is best explained as the result of an orbiting companion. 

The companion to p CrB may be considered to be a "planet", where planets are taken to be stellar 
companions having masses comparable to the mass of the giant planets in our solar system (cf Marcy et al. 
1997). In addition, it is interesting to inquire if there are clues whether this companion was created like 
the planets in our own solar system — that is, by agglomeration of planetesimals in a circularly rotating 
circumstellar disk (e.g. Lissauer 1995), or whether it might have been created in the manner of brown 
dwarfs — that is, by fragmentation within a collapsing molecular cloud, as the low-mass tail of a binary star 
creation. While it has been suggested that a lower limit to the mass of brown dwarfs formed in this way 
is about 10 Mj up due to opacity effects (e.g. Silk 1977), it has recently been argued (Black 1997) that 
uncertainties in the binary formation process allow for brown dwarf companions with much lower mass, 
even overlapping the mass range of giant planets, so by itself the mass of the companion to p CrB does not 
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determine its mode of origin. 

Binary stars, and by extension stars with brown dwarf companions, typically form with significant 
eccentricity. If the orbital period is relatively short, the orbit may be circularized by tidal effects. However, 
main-sequence binaries with periods as long as the 40 day period of p CrB generally do not have circular 
orbits. For example, Latham et al. (1992), in an analysis of old (10-15 Gyr) binaries in the halo, showed 
that only those with periods less than about 19 days were circularized. In addition, an analysis which 
takes into account the effect of tides on the planet itself (Rasio et al. 1996) implies a timescale for orbital 
circularization of the p CrB system of more than 10 12 yr. Thus it appears unlikely that p CrB and its 
companion formed in the manner of binary stars; formation of the companion by agglomeration of dust and 
planetesimals in a dissipative circumstellar disk around p CrB seems more plausible. 

The companion to 47 UMa (Butler & Marcy 1996) is similar to the companion to p CrB, in having 
both near-zero eccentricity and an orbital period long enough that tidal circularization should not have 
occured over the lifetime of the system. However, the companion to p CrB orbits 9 times closer to its parent 
star than does the companion to 47 UMa (at 2.1 AU), and 23 times closer than does Jupiter (at 5.2 AU). 
The companions to both 47 UMa and p CrB lie inside the radius of the ice condensation zone, which has 
been argued (Boss 1995) to be the minimum radius for formation of giant planets. The existence of a giant 
planet, in near-circular orbit as close as 0.23 AU to its parent star but still outside the tidal circularization 
radius, sharpens the question of whether such bodies could be formed in situ, or if not, how they could 
migrate inward and stop in a stable orbit at their present location. 

We are grateful to the Mt. Hopkins observing and support staff, especially Wayne Peters, Bastiaan 
van't Sant, Perry Berlind, and Jim Peters, for help with the observations. We thank Paul Butler and 
Geoffrey Marcy for obtaining very high resolution spectra of our I 2 absorption cell using the McMath-Pierce 
FTS spectrograph. We are indebted to Guillermo Torres and David Latham for useful discussions. We 
thank Gregory Henry for allowing us to quote his photometric results on p CrB in advance of publication. 
We are grateful to the referee, Geoffrey Marcy, for careful reading of the manuscript and helpful suggestions. 
In preparation of this paper, we made use of the Simbad database operated at CDS, Strasbourg, France 
and the NASA Astrophysics Data System. 
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Tabic 1. Orbital Parameters for p CrB 



p 


39.645 ± 0.088 days 


K x 


67.4 ±2.2 ins" 1 


e 


0.028 ±0.040 




210 ± 74 deg 


^pcriastron 


2, 450, 413.7 ±8.2 HJD 


ai sinz 


(36.75 ±0.92) x 10 6 m 


fi(m) 


(1.258 ±0.093) x 10~ 9 M 


^transit 


2, 450, 559.37 ±0.54 HJD 


N 


38 


rms(0 - C) 


9.2 iris" 1 



